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ABSTRACT 

A  simplified  method  for  computing  the  propeller  inflow  velocity  on 
ships  with  propellers  mounted  on  the  side  or  underneath,  such  as  destroyers, 
is  developed.  The  method  utilizes  the  superposition  of  quasi -two-dimensional 
wake  calculations  of  struts  and  shafts  and  the  accurate  calculation  of  the 
nominal  3-D  boundary  layer  and  potential  flow  that  occurs  on  the  bare  hull. 

It  is  possible  to  calculate  the  nominal  flow  for  these  cases  because  the 
propellers  are  mounted  forward  of  the  stern  region.  There  the  thin  boundary 
layer  calculation  methods  are  valid. 
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Section  1 
I  PRODUCTION 


An  analytical  method  for  predicting  the  flow  around  and  in  the  wake 
of  hull  appendages  is  desirable  due  to  the  scaling  difficulties  inherent  in 
appended  model  testing.  The  flow  around  the  appendages  is  Reynolds  number 
(R^)  dependent  and  is  not  correctly  extrapolated  to  full-scale  by  Froude  number 
(F^).  An  analytical  solution  for  the  viscous  flow  around  an  appended  hull  is 

difficult  due  to  flow  separation  around  the  appendages  at  R  of  interest 

a  9  ^ 

(10  <  R  <  10  )  and  the  fuct  that  parts  of  the  appendages  are  embedded  in  the 

"  'n  ^ 

hull  boundary  layer.  It  is  not  yet  feasible  to  solve  the  Navi er-Stokes  equa¬ 
tions  for  flow  around  appendages  on  the  hull  so  that  approximate  techniques 
should  be  developed.  For  certain  hull  appendage  configurations,  such  as 
destroyer  hull  propeller  shafts  and  shaft  support  struts,  approximate  calcula¬ 
tion  methods  can  be  worked  out.  Such  methods  that  can  predict  fairly  accurate 
values  of  the  propeller  inflow  would  be  extremely  valuable  for  design  purposes. 
The  fact  that  on  destroyer  hulls  the  propellers  are  not  completely  immersed  in 
the  hull  boundary  layer  and  that  hul 1 -appendage  interference  effects  are  not 
in  the  way  of  the  propeller  are  reasons  why  simple  approximate  models  for  the 
propeller  inflow  may  be  feasible.  The  main  problem  to  solve  in  this  case  is 
obtaining  an  accurate  description  of  the  nominal  flow  over  the  hull  (including 
wave  effects)  with  the  appendages  removed.  This  nominal  flow  is  needed  to 
trace  the  path  of  the  appendage  wakes  which  are  largely  two-dimensional  and 
can  be  fairly  accurately  calculated.  For  destroyer  hull  forms  it  is  not  neces¬ 
sary  to  deal  with  complicated  stern  shapes,  such  as  cruiser  sterns,  because  the 
propellers  are  well  forward  of  the  stern.  It  is  only  necessary  to  calculate  the 
nominal  flow  over  the  hull  forward  of  the  stern  precisely  in  the  region  where 
current  3-D  boundary  layer  and  potential  flow  methods  work  well. 

For  application  to  ships  whose  propellers  are  mounted  along-side  of 
the  hull  and  ahead  of  the  stern,  this  study  develops  a  simplified  method  for 
computing  the  flow  into  a  propeller.  This  flow  is  only  partially  disturbed  by 
the  hull  boundary  layer  but  can  be  substantially  disturbed  by  the  appendage 
wakes . 
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Section  2 

TECHNICAL  APPROACH 


Recently  SAIC  has  been  engaged  in  making  bare  hull  viscous  flow 
calculations  for  the  DDG  51  destroyer  hull  form  (von  Kerczek,  et  al . ,  1983, 
and  Stern  and  von  Kerczek,  1983).  Figure  1  shows  the  DDG  51  stations  in  way 
of  the  propeller  and  the  shaft  support  struts.  Also  shown  is  the  propeller 
disc  and  the  calculated  nominal  boundary  thickness  at  the  propeller  station 
for  both  model  and  full-scale  and  for  the  30-knot  full-scale  speed  condition. 

It  can  be  seen  in  Figure  1  that  for  model-scale  most  of  the  propeller  disc  is 
outside  the  bare  hull  boundary  layer  and  that  for  full-scale  all  of  the  pro¬ 
peller  disc  is  outside  the  bare  hull  boundary  layer.  It  is  also  evident  that 
the  wake  from  the  propeller  shaft  and  shaft  support  struts  flows  directly  into 
the  propeller  disc.  The  experimental  nominal  wake  results  for  the  DDG  51 
(Dawson,  1982)  confirm  that  most  of  the  wake  deficit  is  due  to  wake  from  the 
appendages.  It  is  also  indicated  by  Dav/son  that  the  wake  velocities  are 
about  10  higher  than  the  model  speed  near  the  lower  part  of  the  propeller 
disc.  This  is  probably  due  to  free  surface  effects  and  implies  the  need  for 
including  such  effects  in  the  bare  hull  potential  flow  calculation.  Also  note 
that  the  strut/hull  intersection  is  about  one  propeller  radius  above  the 
propel ler  disc. 

Based  on  the  above  observations  the  following  model  is  proposed.  It 
is  assui:ied  that  the  primary  effects  of  the  interaction  between  the  bare  hull 
and  the  appendage  boundary  layers  is  that  the  appendage  boundary  layer/wake 
development  is  driven  not  just  by  the  external  potential  flow  but  by  the  local 
velocity.  This  local  velocity  may  be  the  hull  boundary  layer  for  those  parts 
of  ttie  appendage  located  close  to  the  hull  or  the  external  potential  flow  for 
those  parts  of  the  appendage  well  outside  the  boundary  layer.  A  further  ap¬ 
pro/,  i'-uit  ion  is  to  superpose  results  from  two-dimensional  strut  and  shaft 
sectif.ri  wale  calculations  at  different  propeller  radii.  This  is  justified  by 
the  low  value  of  local  flow  sweep  angle  on  the  struts  indicated  by  the  bare 
hull  viscous  flow  solution.  Unfortunately,  funding  constraints  did  not  permit 
free  surface  effects  to  be  included.  However,  future  effort  should  be  directed 
towards  including  free  surface  effects. 


2-1 


Table  I 


WAKE  DEFECT  PREDICTIONS  DUE  TO  STRUTS 


CONDITION 

r/R 

WAKE  WIDTH  AT 
PROPELLER  PLANE 

AXIAL  VELOCITY  RATIO 

AT  WAKE  CENTERLINE 

Full  Scale 

0.63o 

12.3° 

0.82 

Full  Scale 

0.781 

o 

o 

o 

0.82 

Full  Scale 

0.963 

8.1° 

0.82 

Model  Scale 

0.633 

14.3° 

0.73 

Model  Scale 

0.781 

11.6° 

0.73 

Model  Scale 

0.963 

9.4° 

0.73 

0  0.02  0.04  0.06  0.08  0.10  0.12  0.14 


oc  — 


10  — 


12  — 

I 

i 

i  SHIP  q 

Figure  4;  ATHFflA  stations  in  vjay  of  the  propeller  and  the  shaft  support 
struts . 
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The  strut  and  shaft  wake  analyses  presented  in  this  report  were 
compared  to  the  Reed  and  Day  (1979)  data  for  the  longitudinal  velocity  in  the 
propeller  plane.  Comparisons  are  presented  at  propeller  radii  of  r/R  =  0.633, 
0.781,  0.963.  No  attempt  was  made  to  compare  at  the  innermost  radius,  r/R  = 
0.456,  because  of  the  complicated  flow  pattern  and  large  scatter  in  the  data. 

In  order  to  apply  the  wake  analyses,  the  bare  hull  potential -flow  and  boundary- 
layer  solutions  are  required.  The  SAIC  potential-flow  and  boundary-layer  codes 
(von  Kerczek,  1982)  were  run  at  model  and  full-scale  conditions.  The  small- 
crossflow  option  was  used  for  the  boundary-layer  solution.  Calculated  nominal 
boundary-layer  thicknesses  at  the  propeller  station  are  plotted  in  Figure  4. 

As  can  be  seen,  most  of  the  propeller  disk  lies  outside  of  the  boundary  layer. 

In  Figure  5,  comparisons  are  shown  for  the  computed  and  measured  velocity 
profiles  at  x/L  =  0.9.  The  present  velocity  profile  predictions  are  similar 
to  those  given  in  Reed  and  Day  (1979).  That  is,  the  full-scale  profile  is  in 
good  agreement  but  the  predicted  model  profile  is  too  full. 

The  struts  on  the  Athena  were  Navy  EPH  sections  with  a  chord  to 
thickness  ratio  of  six.  For  the  results  of  this  report,  the  struts  were 
modeled  by  an  NACA  0021  airfoil  section.  This  was  for  convenience  because  the 
geometry  and  potential  flow  were  available.  Since  both  sections  are  similar 
in  shape,  it  is  felt  that  this  approximation  will  have  a  negligible  effect  on 
the  predicted  wake  defects.  For  initialization  of  the  wake  calculations.  Head's 
boundary-layer  equations  were  integrated  over  the  NACA  0021  airfoil  section  for 
model  and  full-scale  conditions.  In  both  cases,  the  airfoil  sections  were 
assumed  fully  turbulent.  In  reality,  there  was  probably  considerable  laminar 
flow  (perhaps  one-half  the  airfoil)  on  the  model  section.  For  full-scale  and 
model-scale  conditions,  the  velocity  defect  due  to  the  strut  was  overpredicted, 
as  shown  in  Table  I.  Due  to  the  scaling,  wake  defects  appear  as  spikes  when 
plotted  as  in  Figure  6.  Similar  results  were  obtained  by  Reed  and  Day  (1979). 

A  predicted  wake  profile  at  the  propeller  plane  is  shown  in  Figure  7.  Most  of 
the  velocity  defect  in  the  propeller  plane  is  due  to  the  shaft.  For  the  shaft 
inclined  at  10°  to  the  flow,  the  stream  sees  elliptical  sections  of  c/t  '  6. 
Reynolds  numbers  based  on  shaft  diameter  are 
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Section  4 

SAMPLE  CALCULATION 


Comparisons  were  made  to  the  data  and  computations  of  Reed  and  Day 
(1979).  In  their  experiment,  they  measured  the  wake  and  boundary  layer  on  a 
full-scale  twin-screw  displacement  ship  -  The  Athena.  The  full-scale  wake 
measurements  were  taken  at  four  radii  in  the  propeller  plane  and  the  full- 
scale  boundary-layer  profiles  were  obtained  at  three  longitudinal  locations 
with  and  without  the  propeller  operating.  The  model-scale  nominal  wake  was 
determined  in  a  towing  tank  using  pitot  tubes  while  the  model-scale  boundary- 
layer  measurements  were  made  in  a  wind  tunnel  using  hot  wire  anemometers. 

For  their  calculations,  Reed  and  Day  (1979)  decomposed  the  velocity 
in  the  propeller  disk  as  follows: 

Velocity  =  uniform  stream 

+  perturbation  due  to  hnll  shape 

+  perturbation  due  to  boundary-layer  displacement  effect 
+  viscous  wake  of  struts 
+  viscous  wake  of  shafting. 

They  found  a  negligible  effect  of  the  hull  potential  flow  on  the  propeller 
plane  velocities.  Since  the  boundary-layer  displacement  thickness  increased  the 
hull  thickness  by  less  than  one  percent  of  beam,  the  boundary-layer  displace¬ 
ment  effect  was  neglected.  The  velocity  defect  caused  by  the  struts  was  pre¬ 
dicted  using  an  empirical  far-wake  formulation.  Their  calculated  velocity  defects 
due  to  struts  significantly  overpredicted  the  velocity  defects  that  were 
observed  on  either  model  or  full  scale.  Reed  and  Day  (1979)  estimated  the  wake 
behind  the  shaft  by  a  relationship  between  the  base  pressure  of  an  inclined 
circular  cylinder  to  the  velocity  defect  in  the  wake.  They  predicted  the 
shaft  wake  to  be  less  than  two  percent  of  model  speed  and  one  percent  of  ship 
speed.  These  were  significantly  less  than  measured  values. 
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Section  3 

PROCEDURE  FOR  OBTAINING  PROPELLER  PLANE  VELOCITIES 


The  following  procedure  is  recommended  for  calculating  velocities 
in  the  propeller  plane,  including  appendage  effects. 

1.  Calculate  the  potential  flow  around  the  ship  without  appendages. 

2.  Calculate  the  3-D  boundary  layer  over  the  ship  without  appendages. 

3.  From  a  particular  location  on  the  appendage,  trace  the  stream¬ 
line  to  the  propeller  plane. 

4.  For  a  particular  location  on  the  appendage,  find  u/u^  from  the 
ship  boundary  layer  solution. 


5.  Use  this  value  of  u/u^  to  drive  the  appendage  wake  calculations. 

6.  At  the  propeller  plane 


“  =  “e  <r-r’«ke 

e  scale 


where  the  potential  flow  velocity  u 


is  taken  at  the  propeller. 
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that  on  most  ships,  except  near  the  propeller  hub,  the  propeller  plane  can  be 

considered  to  be  in  the  far  wake  of  the  shaft.  It  is  known  that  in  the  2-D 

1  /2 

turbulent  far  wake  of  cylinders,  the  wake  width  b  n  x  and  the  wake  center- 

-1/2 

line  velocity  u^j^  ax  .  Similarity  solutions  can  be  obtained,  e.g., 
Schlichting  (1968).  Constants  based  on  Schlichting's  (1968)  measurements  for 
wakes  behind  circular  cylinders  of  diameter  d,  give 


where 


b 

W 


0.569  (xCj^d)^/^ 

0.976  (A)’^''^ 
D 


3/2l  2 

1  -  (y/b) 


(22) 


W  =  u  -  u , 
e 


=  — p —  ,  A  =  frontal  area  . 

^  PUg  A 

From  Schlichting's  measurements  these  formulas  are  accurate  for 


50 


Since  the  shaft  is  inclined  at  some  angle  to  the  flow  (typically 
about  10°),  the  drag  coefficient  C„  should  be  for  elliptical  sections.  Hoerner 
(1969)  presents  equations  for  drag  coefficients  of  elliptical  sections.  They  are 


Cp  =  0.015  (1+c/t)  +  1.1  (t/c) 

Cp  =  0.005  (4+2(c/t)  +  120  (t/c)^) 


for  uCg/v  <  1x10^ 
for  uc  /v  >  1x10^ 


(23) 


where  c/t  is  the  chord  to  thickness  ratio.  For  wake  calculations  presented  in 
this  report,  the  virtual  origin  was  taken  at  the  shaft  centerline. 


c 


ave 


I 

L 


(17) 


For  a  flat  plate,  the  Blasius  turbulent  skin-friction  formula  yields 


c  =  1.25  c^(L)  (18) 

‘ave 

where  the  local  value  of  c^  has  been  evaluated  at  the  plat  length  L.  If 


u* 

ave 


ave 


one  obtains 


u*  =  u*  7(1.25)^^^ 
ave 


(19) 


(20) 


for  a  flat  plate.  This  relationship  for  u*,  equation  (20),  is  substituted 

into  Alber's  equations  (13)  and  (15).  It  is  assumed  that  the  use  of  u*  in 
^  ave 

Alber's  (1980)  equations  is  a  good  representation  for  ^Qi_for  airfoil  wakes, 
although  the  authors  were  unable  to  obtain  data  for  comparisons. 


Once  the  wake  centerline  velocity  and  wake  width  have  been  calculated, 
a  wake  velocity  profile  can  be  obtained.  Ramaprian  et  al.  (1982)  found  that 
the  velocity  defect  profile  in  far-wake  coordinates 


EXP 


n  =  y/b 


(21) 


accurately  represents  wake  data,  even  for  x/y  =  25.  This  formulation  was  con¬ 
sidered  Sufficient  for  the  present  applications.  In  equation  (21),  W  =  u^-u 
is  the  velocity  defect,  is  the  maximum  velocity  defect,  and  b  is 
the  wake  hal f  wi dth . 


The  flow  over  an  inclined  shaft  is  fully  separated  and  thus  is  com¬ 
putationally  difficult.  Such  a  calculation  was  beyond  the  scope  of  this  task. 
Instead  the  following  approximate  procedure  was  employed.  First,  it  is  noted 
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where  is  the  wake  centerline  velocity  and  u*  is  the  friction  velocity  at 
the  plate  trailing  edge.  Ramaprian  et  al.  (1982)  show  that  the  equation  sug¬ 
gested  by  Alber  (1980)  is  more  accurate  when  compared  to  the  data  sets  of 
Ramaprian  et  al.  (1982),  Pot  (1979),  Andreopoulos  (1980)  and  Chevrey  and 
Kovasznay  (1969).  Alber  (1980)  obtained  an  analytical  solution  using  an 
inner  and  outer  layer  coordinate  expansion  technique.  It  appears  accurate  to 
X*  greater  than  2x10  .  Alber's  (1980)  equation  is 


^  =  j  [>ng(x*)  -  y|  +  B  (13) 

where  K  =  0.41,  y  =  0.5772157  (Euler  constant),  B  =  5.0.  The  function  g  is 
defined  by 


g(x*) 


In  g(x*) 


(14) 


In  the  central  portion  of  the  wake  immediately  downstream  of  the  trailing 
edge,  Alber  assumed  the  flow  to  be  laminar  and  recovered  the  Goldstein  solution 


=  1.61  .  (15) 

u* 

This  solution  is  valid  until  the  inner  wake  grows  to  equal  the  upstream  sub¬ 
layer  thickness.  The  streamwise  extent  of  this  region  is  on  the  order  of  ten 
sublayer  thicknesses.  Equations  (13)  and  (15)  match  at  approximately  x*  =  300, 
so  this  will  be  the  crossover  criterion  used.  In  its  present  form,  the  u^j^ 
calculation  is  based  on  u*  determined  from  the  local  c^  evaluated  at  the 
plate  trailing  edge  from 


=  (B/Cf)'/'  .  (16) 

For  an  airfoil  this  is  a  problem  since  c^,  and  hence  u*,  tends  to  zero  very 
rapidly  near  the  trailing  edge.  To  avoid  this  problem,  Alber's  (1980)  cor¬ 
relation  is  evaluated  with  u*  based  on  c^  where 

a  ve 
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Figure  3:  Wake  shape  factor  comparison. 


known  from  measurements.  Adjustment  of  the  wake  entrainment  rate  follows 
the  work  of  Green  (et  al . ,  1973).  From  Green's  analysis,  just  downstream 
of  the  trailing  edge 


6  =  -0.06  (H-1)^  , 

dx 


while  in  the  far  wake 

e  ~  =  -0.242  (H-1)^  .  (9) 

dx 

Here  6  is  the  half-wake  momentum  thickness.  Based  on  the  findings  of 
Ramaprian  et  al.  (1982),  flat  plate  far  wake  conditions  are  obtained  when  x/0 
is  greater  than  350.  It  was  decided  to  adjust  the  wake  entrainment  by  cal¬ 
culating  it  from  an  equation  of  the  form 

0  ■—  =  C  (H-l)^  (10) 

The  constant  C  is  linearly  interpolated  from  its  trailing  edge  (t.e.)  value 
(-0.06)  to  its  far  wake  value  (-0.242),  by 

-0.06  (350  -  X  /9)  -  0.242  (x  /B) 

r  =  _ ° ° _  (11) 

^  350  ^  ‘ 

where  x  =  x  -  x^  .  Comparison  of  the  calculated  H  to  the  flat  plate  wake 
0  t.e. 

data  of  Chevrey  and  Kovasznay  (1969)  is  shown  in  Figure  3.  Agreement  is  seen 
to  be  excellent.  Integration  of  equations  (1),  (2)  and  (10)  (with  c^=0)  along 
with  equation  (7)  provides  the  wake  width.  Additional  calculations  are  needed 
to  obtain  the  wake  velocity  profile  and  centerline  velocity. 

Both  Andreopoulos  and  Bradshaw  (1980)  and  Alber  (1980)  have  suggested 
wake  centerline  velocity  correlations  of  the  form 


Ucl/u*  = 


2-5 


The  boundary- 1 ayer  entrainment  equation  is 


dx 


F(G)  -  S  ^ 
u  dx 


(2) 


where 


G  =  1.535  (H-0.7)  +  3.3 


(3) 


and 


F(G)  =  0.0306  (G-3.0) 


-0.653 


The  skin- friction  coefficient  was  modeled  by  the  Ludwieg-Ti llmann 


c  =  0.246  X 

T  U 


but  any  other  =  c^(Rey,H)  could  have  been  used.  Since 


G  =  (6-6*)/e,  q  =  G9 


one  obtains 


6  =  q  +  OH  . 


(4) 

formula 

(5) 


(6) 


(7) 


On  many  airfoils  the  flow  separates  near  the  trailing  edge.  Rather 
than  attempt  to  calculate  through  separation,  the  wake  calculations  are  started 
at  the  location  of  separation.  The  turbulent  wake  from  the  airfoil  trailing 
edge  to  the  propeller  plane  is  calculated  using  basically  the  same  method  as 
for  the  turbulent  boundary  layer  with  the  following  additions.  The  wake 
centerline  is  assumed  to  follow  the  bare  hull  potential  flow  streamline  from 
the  airfoil  trailing  edge  to  the  propeller  plane.  In  the  wake  the  skin-friction 
is  zero  and  the  wake  entrainment  rate  is  adjusted  to  approach  its  far  wake  value 


The  calculation  of  the  boundary  layer/wake  of  a  strut  (or  shaft) 
section  is  a  four-step  procedure  (see  Figure  2), 

(1)  Calculate  laminar  boundary  layer  from  the  strut  leading  edge 
stagnation  point  to  the  point  of  transition  to  turbulent  flow, 

tr 

(2)  Determine 

(3)  Calculate  turbulent  boundary  layer  from  to  the  strut 
trailing  edge,  x^^ . 

(4)  Calculate  turbulent  wake  from  x^^  to  propeller  plane. 


Figure  2;  Sketch  of  strut-section  boundary  layer/wake. 

For  the  present  calculations  the  flow  over  the  appendages  has  been 
assumed  all  turbulent.  Two  types  of  appendages  are  considered:  struts  and 
shafts . 


The  turbulent  boundary  layer  over  the  struts  is  calculated  using  Head's 
(1958)  momentum  integral  entrainment  method.  For  a  2-D  flow,  the  momentum  inte¬ 
gral  equation  is 
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Fiaure  6;  Uake  streamwise  velocity  ratio  for  ATHENA  -  full  scale  strut 
effect. 


Figure  7:  Typical  wake  profile  prediction. 


R.  =  ud/v  =  1.6  X  10  full  scale 

4 

=  6.7  X  10  model  scale. 

Then,  formulas  (23)  give 

Cp  =  0.1  full  scale 

Cp  =  0.3  model  scale. 

The  wake  half  width  and  velocity  defect  were  calculated  from  equation  (22). 
These  results  are  plotted  in  Figures  8-10  for  full  scale  and  Figures  11-13 
for  model  scale.  The  wake  widths  seem  correct  for  both  model  and  full  scale, 
but  the  wake  defect  for  the  model  scale  is  too  large.  Of  course,  all  the  data 
lie  above  the  predictions  because  of  free  surface  effects. 


4-7 


O  FULL  SCALE  DATA  (REED  AND  DAY) 
X  r.  y.  MODEL  SCALE  DATA 


PRESENT  PREDICTION  (FULL  SCALE) 


3  3  20  *0  SO  3C  100  120  MO  ISO  1 8C  200  220  240  260  290  300  320  3*0  360  390 

ANCLE  IN  DECREES 


Figure  8:  Wake  streamwise  velocity  ratio  for  ATHENA  -  full  scale 
shaft  effect  at  0.633  radius. 
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Figure  9;  Wake  streamwise  velocity  ratio  for  ATHENA  -  full  scale 
shaft  effect  at  0.781  radius. 
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Figure  12:  Wake  streamwise  velocity  raito  for  ATHENA  -  model 
scale  shaft  effect  at  0.781  radius. 
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Figure  13:  Hake  streamwise  velocity  ratio  for  ATHENA  -  model 
scale  shaft  effect  at  0.963  radius. 
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Section  5 
CONCLUSIONS 


Based  on  the  available  data  it  is  not  possible  to  assess  fully 
the  effectiveness  of  the  appendage  wake  calculation  method  presented  here. 

The  comparisons  between  the  calculated  and  measured  wakes  do  show  fair  agree¬ 
ment.  However,  the  measurements  do  seem  sufficiently  erratic,  especially  the 
full  scale  ones,  that  better  agreement  can  hardly  be  expected.  Since  the 
appendage  affect  on  propeller  inflow  is  not  large,  the  fair  agreement  obtained 
by  the  calculation  method  may  be  adequate,  at  least  for  the  purpose  of  powering 
predictions.  However,  for  the  purpose  of  predicting  vibration  and  noise 
characteristics  of  the  propeller,  where  the  harmonic  distribution  of  wake 
defect  is  important,  more  accurate  procedures  may  be  desirable. 

What  does  seem  clear  in  the  experiments  is  that  the  propeller  inflow 
is  mainly  affected  by  the  potential  flow  and  geometry  of  the  hull  in  the 
vicinity  of  the  propeller  disks.  It  seems  important  to  include  wave  effects 
since  the  local  velocity  is  substantially  affected  by  waves.  Unfortunately 
funding  for  this  project  did  not  include  means  to  calculate  the  full  free 
surface  wave  effects  on  propeller  inflow. 

It  seems  worthwhile  to  obtain  more  experimental  data  on  appendage 
wake  effects  to  further  validate  and  improve  the  tool  developed  in  this  effort. 
Once  the  present  tool  is  proved  to  be  uniformly  reliable  to  at  least  the 
accuracy  indicated  here  then  it  should  prove  useful  to  designers  for  evalua¬ 
tion  purposes. 
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